NSI=size(M0,1);
[XC,YC]=meshgrid((1:NSI)-NSI/2-1);
WRC=sqrt(XC.^2+YC.^2)/(pseq.fovcm);
MRC=WRC<(1.128*2.5);
R2g=abs(fftshift(ifft2(fftshift(MRC.*fftshift(fft2(fftshift(R2e)))))));
M0g=abs(fftshift(ifft2(fftshift(MRC.*fftshift(fft2(fftshift(M0)))))));